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I. INTRODUCTION 

At high energy or, equivalently small x, the interac- 
tions of hadrons are expected to be dominated by non- 
linear strong color fields. A convenient effective theory 
approach to studying these color fields is provided by the 
Color Glass Condensate (for reviews see e.g. [UH]). This 
effective theory is based on the division of the QCD de- 
grees of freedom into a static (in light cone time) large x 
part, which is described as a color source, and the small x 
gluons as a color field radiated by these sources. The ra- 
pidity scale separating the two is, of course, arbitrary, 
and should be chosen to correspond to the rapidity scale 
probed in the particular physical process being looked 
at. From the requirement that physical cross sections 
must not depend on this arbitrary scale one can derive 
renormalization group equations describing the rapidity 
dependence of the probability distribution of the color 
sources. The general, leading In a;, equation for the prob- 
ability distribution of different color charge configura- 
tions is known by the acronym JIMWLK (see Ref. [5] 
and references therein) . The mean field approximation of 
JIMWLK, known as the Balitsky-Kovchegov [4-6] equa- 
tion, is widely used in phenomenological applications. 
The BK equation describes the rapidity dependence of 
the unintegrated gluon distribution and can loosely be 
thought of as describing the rapidity dependence of the 
mean number of gluons in stead of the whole probability 
distribution encoded in the JIMWLK equation. 

The BK equation is usually enough when studying sin- 
gle inclusive cross sections in a dilute-dense collision such 
as DIS or proton-nucleus collisions at forward rapidities. 
Understanding the initial state of a heavy-ion collision, 
on the other hand, necessitates the, eventually numer- 
ical, solution of a much more complicated interacting 
system [7J [S] . Therefore these very inclusive quantities 
do not provide direct experimental access to all of the 
physics described by the JIMWLK equation (for a discus- 
sion of potential observables in diffractive DIS that go be- 
yond the mean- field BK approximation see e.g. Ref. [§])• 



It has recently become more evident that multiparticle 
correlations could provide precisely such an observable. 

Quite generically multiparton correlations are more 
sensitive to the detailed dynamics of the colliding ob- 
jects than single particle distributions. An attractive ob- 
servable, measurable at both RHIC and the LHC, are 
dihadron correlations in the forward rapidity region in 
proton-nucleus or deuteron-nucleus collisions. In this 
kincmatical regime the dilute and relatively well under- 
stood large x part of the proton provides with a cali- 
brated probe of the nonlinear small x gluon fields in the 
target nucleus. By varying the rapidity one can study 
different values of x and the transverse momentum de- 
pendence of the correlation should be directly sensitive 
to the dominant intrinsic transverse momentum scale in 
the target, the saturation scale Q s . 

Indeed measurements of correlations between two for- 
ward dihadrons measured in dAu collisions at RHIC jTOl 
ITT] seem to show indications of "initial state" or "cold nu- 
clear matter" effects that are significantly stronger than 
in pp collisions or at central rapidities [LWH] . The up- 
coming LHC proton-lead collisions will provide more op- 
portunities to study these phenomena in a wider kine- 
matical range. These observations have provided an im- 
petus for renewed interest in the gluonic correlations in- 
cluded in the JIMWLK evolution [T5HH]. In particular 
it was argued HMS] that the result of a full JIMWLK 
evolution, also at finite N c , can quite accurately be cap- 
tured by the so called Gaussian approximation, relat- 
ing higher point Wilson line correlators to the two-point 
function. Thus in the Gaussian approximation one can 
construct observable cross sections using the solution of 
the BK equation alone, which is much more convenient in 
practice than solving the full JIMWLK equation. These 
recent theoretial developments were not fully reflected 
in the pioneering calculations of dihadron correlations 
in [221 123] (see also the more recent work [24 , 25j where a 
fc^-factorized approximation is derived in a certain kine- 
matical limit). The main purpose of this paper is to 
implement the Gaussian approximation, which so far has 
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only been tested for particular coordinate space configu- 
rations, in a full calculation of the dihadron correlation 
We shall first, in Sec. [TTJ discuss the description of the 
single inclusive hadron spectrum as a baseline that should 
be consistently described by the same calculation and 
used to constrain the initial conditions of evolution in 
rapidity. We then, in Sec. |III[ recall the expression of the 
dihadron cross section derived in Ref. using the light 



cone perturbation theory formalism. In Sec. IV we show 
that the dihadron cross section as derived in Ref. [22] con- 
tains a logarithmically infrared divergent part, which has 
been overlooked in the previous literature. We identify 
this as a double parton scattering contribution describ- 
ing the independent scattering of two partons already 
present in the dilute projectile. The double parton scat- 
tering contribution must be consistently subtracted from 
the correlated cross section and absorbed into a separate 
contribution that depends on additional nonperturbative 
information about the projectile in the form of a dou- 
ble parton distribution. We then discuss in Sec. [V] the 
expressions for the Wilson line correlators used in this 
work. In Sec. |VI| we present our results for the dihadron 
cross section. We show that using the full expression for 
the dihadron cross section enhances the away-side peak 
by a large factor compared to approximations used in 
the previous literature. The qualitative features of the 
results in the previous literature remain, most promi- 
nently a strong dependence of the correlation on px/Q s , 
which leads to a large difference between deuteron-gold 
and proton-proton collisions. Technical details on the 
DPS limit are relegated to Appendix [5] and on the im- 
pact parameter dependence to Appendix [B] 



II. BASELINE DESCRIPTION OF THE 
TARGET: BK EVOLUTION AND SINGLE 
INCLUSIVE SPECTRA 

We shall here use the dipole cross sections obtained 
from solving numerically the BK evolution equation us- 
ing the Balitsky running coupling prescription [28] (see 
the comparison of different running coupling schemes 
in Ref. [29] ). The evolution equations must be supple- 
mented with an initial condition at a starting rapidity (or 
x). Ideally one would like to fully constrain the initial 
conditions of BK/JIMWLK evolution by a comparison 
to small x DIS data. Due to the large amount of precise 
data available from HERA, the expectation value of the 
dipole cross section in a proton is quite well constrained. 
Several ways to extend these parametrization to nuclei 
exist in the literature, and in practice it is useful to also 
compare with single inclusive hadron production data 
in pA and dAu to further constrain the nuclear dipole 
cross sections. We shall here first recall the ("hybrid for- 
malism" ) equations used to calculate the single inclusive 
hadron yields at forward rapidity. We then discuss the 
MV model parametrization of the initial condition used 
in this paper and compare it to the experimental forward 
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FIG. 1: The identified hadron spectra at forward rapidity 
compared to the experimental data from proton-proton colli- 
sions from the STAR and BRAHMS [27] collaborations. 
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FIG. 2: The nuclear modification factor i?dAu using our 
MV dipole cross section parametrization. The experimental 
results shown are PHENIX [11] centrality dependent (solid 
curves, 3 < y < 3.8) and BRAHMS [27] {y = 3.2, upper 
dashed line) and STAR [2B] (y = 4, lower dashed line) mini- 
mum bias results. The calculations are the spectra obtained 
using Eq. S. 



single inclusive spectra. 

The single inclusive yields for quark and gluon produc- 
tion are given by 



dy d 2 q T 

c Y N hA^gX 

dy d 2 q T 



1 



(2tt)= 
1 



xq{x)S{q T ), 
xg(x)S(q T ), 



(1) 
(2) 
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where S(qr) is the Fourier-transform of the fundamental 
representation dipole operator (D) and S(qr) the corre- 
sponding quantity in the adjoint representation. Here 
xq(x) and xg{x) are standard collinear parton distribu- 
tion functions describing the large x partons from the 
dilute projectile, for which we use the CTEQ NLO [55] 
parametrization. To get the single hadron spectrum we 
convolute the parton-level cross section with the DSS [53] 
fragmentation function, resulting in 

dN hA^h'x r dz dN hA ^ tX . 2n , . 

i A2 = y, / ~T a A2i rx D i-+h'{z,IJ> )■ (3) 

dyd 2 q T J z 1 dyd 2 {q T /z) 

Throughout this work the scale fj, at which the parton 
distribution function, fragmentation function and strong 
coupling constant a s are evaluated is chosen to be the 
transverse momentum of the produced hadron in single 
inclusive production and transverse momentum of the 
leading hadron in dihadron production. 

As an initial condition we have used the MV model 
parameters used in Ref. [35] to fit single inclusive hadron 
production data at RHIC. The initial dipole cross section 
is 

with an initial saturation scale Q 2 = 0.2 GeV 2 at xq = 
0.007 for the proton, as in Ref. [55] . A comparison of the 
resulting yields in proton-proton collisions to the forward 
results of the STAR [25J and BRAHMS [27] experiments 
is shown in Fig. [IJ While the description of the data 
is not perfect, we are confident that this is a realistic 
enough baseline for understanding nuclear modifications 
to the dihadron yield. Note that there are no arbitrary 
normalization if -factors used in this work. 

Estimates based on nuclear geometry would sug- 
gest [55] [37] that the initial saturation scale Q 2 should be 
even as large as A 1 / 3 ~ 6 times that of the proton (for 
minimum bias collisions). Based on both experimental 
data and the general expectation that saturation effects 
should vanish for px 3> Q s , it would also seem natural for 
the nuclear modification factor i?dAu to approach 1 for 
high pt at least at midrapidity. Assuming that the trans- 
verse area of a proton also at small x is identified with the 
total inelastic nucleon-nucleon cross section, this would 
happen with Q 2 A ~ N^Q 2 p ps 3.6Q 2 p for a minumum 
bias dAu collision 1 . 



One would generically expect Qg to be proportional to the 
number of overlapping nucleons, a quantity which is in prac- 
tice encoded as the number of binary collisions in a Glauber 
model for a proton-nucleus collision. In this paper we will as- 
sume = N££/2 for simplicity. Note that AT^A is also the 
number of participant nucleons in the target nucleus in a pA 
collision. 



However, fits to identified hadron spectra, especially 
at forward rapidities tend to favor a smaller nuclear sat- 
uration scale. In the calculation of Ref. 35j, for exam- 
ple, a calculation of hadron spectra in deuteron-gold col- 
lisions is made with an initial nuclear saturation scale 
Q 2 A s» 2.3Q 2 p (comparing Q s a and Q sp at x = 0.007, 
taking into account the different Xq used for protons and 
for nuclei). These saturation scales give a good descrip- 
tion of the BRAHMS charged hadron spectra at 77 = 2.2 
and 77 = 3.2. 

For the most forward ir° data the nuclear modification 
ratio i?dAu seems to decrease faster with y in the data 
than predicted by running coupling BK evolution. For 
example, in the calculation of Ref. [35] the STAR very 
forward rj — 4 data requires an additional suppression by 
if-factors 0.4 for pp and 0.3 for dAu and a similar effect 
is seen in [57] [38] . This is problematic for the dihadron 
correlation, for which the STAR forward data comes pre- 
cisely at this same kinematical region. 

In view of these uncertainties, we shall not attempt 
to provide the best possible fit to all the existing sin- 
gle inclusive particle production data in this work. In- 
stead we use as an initial condition for nuclei the same 
parametrization Q as for protons, with a nuclear sat- 
uration scale obtained by straightforward iVbi n scaling. 
This is essentially equivalent to the Monte Carlo rcBK 
model of Ref. [59] , where the saturation scale Q 2 is taken 
to be proportional to the number of overlapping nucle- 
ons in a position in the transverse plane. In particular 
it leads, by construction, to a nuclear modification fac- 
tor i?dAu approaching unity at high transverse momenta, 
when the unintegrated parton distribution is probed at 
the scale of the initial condition at xq. For minimum 
bias collisions this leads to an initial nuclear saturation 
scale Q 2 — 0.72 GeV 2 , and for a 0-20% centrality class 
Q 2 = 1.51 GeV 2 , using the N^ia values from Ref. [IT] . 

We show in Fig. [2] the resulting forward ir° nuclear 
modification ratio i?dAu compared to the PHENIX [TT] 
centrality dependent and BRAHMS [27] and STAR [25] 
minimum bias data. Let us now point the reader to some 
features in this plot. Firstly, none of the nuclear data 
compared to here has been used as an input in the calcu- 
lation; in stead the initial condition for BK evolution has 
been set from a simple and straightforward Nbw scaling 
of the saturation scale. Secondly, the STAR minimum 
bias data exhibits approximately as much suppression as 
the most central PHENIX data. Although the two are 
at slightly different rapidities, it seems quite generically 
impossible to simultaneously describe both data sets in 
a unique parametrization as simple as ours. A more de- 
tailed treatment of minimum bias collisions by averag- 
ing over different centralities in the calculation (in stead 
of a single Q s value characterising minimum bias colli- 
sions) could also be expected to improve ths description. 
Note that the negative hadron nuclear modification ratio 
becomes > 1 at high pt due to the difference between 
the deuteron and proton probes (a relative excess of d 
quarks in the deuteron over the proton leads to an en- 
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hancement in negative particle production independently 
of the target). All in all, given that this is a pure pre- 
diction for the nuclear modifiction of single particle pro- 
duction, without arbitrary K factors or adjusting the pa- 
rameters to the nuclear data, we consider this as an ade- 
quate parametrization for the purposes of understanding 
the basic features of dihadron correlations, which is the 
main purpose of this paper. 

III. DIHADRON CORRELATIONS 

We shall consider the scattering process for forward 
dihadron production as a large x quark with momen- 



tum p + from the probe deuteron or proton, propagating 
eikonally through the target nucleus or proton. It can 
radiate a gluon with momentum k + = zp + and is left 
with a longitudinal momentum q + = (1 — z)p + . In the 
high energy limit the scattering of both the quark and 
the gluon can be described by an eikonal approximation, 
where they pick up a phase given by a Wilson line in the 
color field of the target. The detailed derivation of the 
double inclusive cross section is performed in Ref. [25] 
and results in the following expression for the qA — > qgX 
cross section: 



dfc+d 2 k T dg+d 2 q T asUF ° [P K q > J (2vr) 2 (2tt) 2 (2^) 2 (2^) 2 

x^^K x t ~ b' T )0^(x T - b T ){^ 4 >(b T ,x T ,b^) - 5( 3 )(b T ,x T ,z^) (5) 

-s^(z T y T ,h'T) + s {2 \zT,z' T )}, 

with zt = zx.t + (1 — z)hT and likewise, z' T = zx' T + (1 — z)h' T . Here the target is described by expectation values 
of the Wilson line operators 



S*W(b T ,b^,x T ,x^) 
S^(b T ,x T ,z^) 
S i2 \z T ,z' T ) 



N/-1 



(Tr (V{h T )V\h' T )t d t c ) [[/(x x )f/t(x^)] cd ) 



S—^ (Tr (vHz' T )t c V(h T )t d ) C/ cd (x T )) 
(Tr (V(z T )V^z' T ))). 



(6) 
(7) 
(8) 



The momenta of the produced gluon and quark are k and 
q respectively. Likewise, xr,Xy should be interpreted 
as the transverse position of the gluon, br,b^ of the 
quark after the scattering and zt, z' t of the quark before 
the scattering; in the amplitude and the complex conju- 
gate respectively. The wave function ip describes q — > qg 
splitting in coordinate space, and its expression in the 
massless limit is given in Eq. (13). In our numerical cal- 
culations we set the quark mass m q = 0.14 GeV, but the 
finite quark mass has little effect on the final results. 

The scattering amplitude is a sum of two terms corre- 
sponding to the radiation of the gluon happening before 
and after the interaction with the target. Out of the three 
operators ([tJ) , (jij) in the gluon is radiated before 
the interaction and in after, with corresponding 
to the interference between gluon radiation before and 
after the target. 

The operators Q, (|8| are expressed in terms of 
fundamental and adjoint representation Wilson lines, de- 
noted here by V and U, respectively. In the following we 



denote the dipole and quadrupole operators by 



D(x T -y T ) = -^Tr^(x T )yt( yT ) (g) 



Q(x T ,y T ,u T ,v T ) = i r ftV(x T )V+(y r )V(u r )V^(v r ), 

(10) 

and their expectation values as S = = (D) and 

Q = (Q). In terms of these the higher point correlators in 
Eq. ([H can be written in the form quoted in Refs. [TKirT7] 



SW(b T ,b^,x T ,x^) = / - 1 D(bT,b' T ) 



y , _ -D(x T , x^)Q(b T , b' T , x' T , x T ) ) (11) 
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s^Hb 



T, XT, ZtJ 



Nr. 



N 2 _ £-°( bT ' x t)D(x t , z' t ) 



Nr. 



-D(b T ,z T ). 



(12) 



We will discuss explicitly in Sec. [V] the Gaussian approx- 
imation used in this paper to obtain the Wilson line cor- 
relators in Eqs. (11) and (12). Computing the dihadron 



correlation in LHC kinematics will require taking into ac- 
count also the gluon-initiated channel, because the probe 
x will not be as large. This would require calculating an 
eight-point function of Wilson lines, which we leave for 
future work. 



IV. DOUBLE PARTON SCATTERING 
CONTRIBUTION 



In the massless limit the wave function product ap- 
pearing in Eq. ([5J is given by 



a/3\ 



&,(u T )<£* (u T ) = 



87T 2 Ut ' Uj 



k+ |u T | 2 K 



r(l + (l-z) 2 ), 



br and 



(13) 



section one 



must also in- 



with ut = xt 

In the full dihadron cross 
elude a "double parton scattering" (DPS) contribution 
that corresponds to the independent scattering of two 
partons from the probe. This contribution, whose impor- 
tance for understanding the Atp-independent pedestal of 
the double inclusive cross section, was emphasized in |40j . 
We shall show that Eq. ([5J contains a part of this (DPS) 
contribution, which must be subtracted from it to avoid 
double counting. This is in contrast with Ref. [2H 122] 
where, in stead of the full cross section Eq. (51), one con- 
siders only the "correlation limit" in which the DPS con- 
tribution vanishes. As discussed in Ref. [23] the full cross 
section can differ significantly from the correlation limit 
when the transverse momenta of the produced particles 
are comparable to the saturation scale, which is precisely 
the regime that is interesting for observing the effects of 
large color fields in the target. 

In the limit where the gluon is far away from the quark: 
|b T - b' T \ ~ |x T - x T | - 1/Q S , u T = \b T ~ x T | > 1/Q S 
the Wilson lines of the quark and the gluon are uncorre- 
cted. Therefore the expectation values in factorize 
into a product of an adjoint representation dipole at the 
location of the gluon and a fundamental representation 
one at the location of the quark. We call this kinemati- 
cal regime the "DPS" limit since, as we will argue in the 
following, it corresponds to a quark and a gluon, present 
already in the wavefunction of the incoming dilute projec- 
tile, scattering independently off the target. Factorizing 
the expectation values and using the fact that the expec- 
tation values must be color singlets we can write in 



the DPS limit as 



S< 4 > (br, x T , b' T , x T ) « S^ s (b T , x T , b' T , x T ) 



DPS 



N c 2 - 

N c . 



-Tr ((V(b T )VHW T )t d t c ))([U^ T )U\^ T )r d ) 



v 2 _ 1 \D(b T , b' T ) )( D*{x T , x' T ) - ^ ) , (14) 



where 



5(x T - x T ) = 



N r 



(15) 



can be identified as the two point function in the adjoint 
representation, normalized to S(xt = x T ) = 1. The 
cross terms vanish in the DPS limit. This is easily 
seen by noting that z is finite as it is fixed by the final 
state kinematics and thus all three coordinates br,xr 
and z' T (or b T ,x T and zr) are much more than 1/Q S 
apart. Consequently factorizes into a product of 
three single Wilson line expectation values, which are all 
zero. The Wilson line operator corresponding to gluon 
emission after the target, S^ 2 - 1 , is finite in the DPS limit. 
It does not, as we show in Appendix [XJ give a divergent 
contribution to the double inclusive cross section even for 
massless quarks. 

In the massless case the integral (|5j) diverges logarith- 
mically in the DPS limit, because is nonzero and 
the wave function product has a large power law tail at 
large ut,u t . This logarithmic divergence is shown ex- 
plicitly in Appendix [A] Physically this means that the 
quark emits a very small transverse momentum gluon. 
The quark and gluon subsequently scatter independently 
off the target. A simple lifetime argument shows that 
the emission of the quark happens Ax + ~ z(l — z)p + ux 2 
before the interaction with the target. Thus the contribu- 
tion from ut — > oo corresponds to a splitting happening 
a long time before the interaction with the target. This 
logaritmically divergent contribution must be regulated 
by confinement scale physics in the wavefunction of the 
projectile. It is in fact exactly the kind of contribution 
that is represented by double parton scattering [4T| [42] . 
The confinement scale physics of the correlations in the 
large x projectile is not calculable in our formalism. We 
shall in stead add it as an additional explicit DPS con- 
tribution. To avoid double counting we must subtract 
the equivalent contribution from our Eq. |5j, which will 
render it manifestly IR finite. We do this by replacing 
5( 4 )(b T ,x T ,b T ,x T ) in Eq. (g by 

^(b T ,x T ,b T ,x T ) = SW(b T ,x T ,b T ,x T ) 

1 \ „ f 1 



X T 



b^ 



Aqcd / V Aqcd 

x5 D 4 ^ s (b T ,XT,b T ,x T ) (16) 



with S'ppg given by Eq. (|14|) 
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After convoluting with the valence quark distribution dependent part of the double differential yield becomes 
in the probe deuteron or proton this correlated, A.(p- 



j pTh a h b ^qgX 

sub 

dy q d 2 k T dy g d 2 q T 



a s C F xq{x,Q 2 ){l - z)(l + (1 - zf) 



1 



d 2 x T d 2 x T d 2 b T d 2 b^ 



Si J (2tt) 2 (2tt) 2 (2tt) 2 (2tt) 5 



(x' T -x T ) pt q T -(b' T -b T ) H 2 (X T - b T ) ■ (x' T - b' T ) 

1 ' ' b^) 2 



{S s ( u >(b T ,x T ,b^) 



(x T - b T ) 2 (x T - u Tl 
S {3) (b T , x T , z' T ) - 5 (3) (zt, x T , b^) 



5( 2 )(z T ,z^)}. (17) 



Here we have moved from a cross section to a differential 
yield by dividing with the transverse area of the target 
S± . Because the integrand is translationally invariant as 
a whole, one out of the four integrals over the transverse 
plane in Eq. ( 17 1 is performed trivially, canceling the fac- 
tor 1 / S± (see also the discussion on the impact parameter 
dependence in Appendix [B]). We assume the Wilson line 
correlators to be independent of the overall impact pa- 
rameter; thus the dependence on centrality only comes 
through the Wilson line correlator, which should corre- 
spond to an average one for the desired centrality class. 

Now that the DPS contribution has been subtracted 
from the dihadron correlation, we must include it sepa- 
rately to get an estimate of the Ac/?-independent back- 
ground. The logarithmic divergence in the DPS limit is 
physically regulated by confinement scale physics, and 
must be absorbed into a new nonperturbative input de- 
scribing the probe, namely the double parton distribution 
function (DPDF). For deuteron-gold collisions we divide 
the deuteron DPDF into two separate parts. The first 
one corresponds to taking two partons from the same 
nucleon in the deuteron, described by a single-nucleon 
double parton distribution Dij{xi,Xj,Q 2 ). We imple- 



ment the kinematical constraint Xi 
Ref. [ID] by taking 



Xj < 1 following 



D\i (xi ,Xj,Q ) — 



fi{xi)fj \ T 



+ fi 



1 — i, 



fj( x j) 



(18) 



where i and j denote the parton species (g,u,d). For 
the deuteron one must sum over the proton and neu- 
tron: Dfj ( Xi ,Xj,Q 2 ) = D*. {xi ,x ]7 Q 2 ) + D% [ Xi , Xj , Q 2 ) 
The second contribution involves taking one parton from 
the neutron and the other one from the proton, which is 
not bound by the same kinematical constraint. In this 
case one must sum; not symmetrize; over the different 
combinations of different partons from different nucleons: 



D$( Xi , x j: Q 2 ) = [fi{ Xi )% ( Xj ) + ff (x t ) /?(*,•)] 



The double parton scattering contribution is the sum of 
these two 



diV, 



h a h b ^ij X 
DPS 



1 



dy 9 d 2 k T d % d 2 qT (2ir) 4 



+ D^\ Xl ,x 3 ,Q 2 ) )S i (k Ti )S j (k Tj ), (20) 



where the two point functions ^(k^i)? SjfeTj) are taken 
in the fundamental or adjoint representation according 
to the parton species 

The total quark-gluon production cross section is now 
the sum of pi and (fl7|): 



dN h a h b ^ijX = dN hah b ^ijX + dN hah b ^ijX_ ( 21 ) 

After convolution with fragmentation functions the 
cross section for double inclusive production of hadrons 
1 and 2 becomes 

dN h « h »^ hih2X r dzi dz 2 

dyi d 2 p T i dy 2 d 2 p T 2 J z\ z\ 

d ]yh a h b ^ijX 

dyi d 2 (pTi/^i) dy 2 d 2 (p T 2/z 2 



-D hul (z 1 )D h2j (z 2 ) (22) 



where i and j are summed over gluons and light quark 
flavors; for the combination i,j — q,g and i, j = g,q the 
parton level cross section includes the correlated cross 
section d ]yhah b ^>-ijX £ Qr ^ e ^j^g om y jjpg 

contribution. This is our final expression that will be 
compared to the experimental yield. We shall now turn 
to the calculation of the different Wilson line correlators 
appearing m diV sub J . 



V. MULTI-POINT CORRELATORS OF WILSON 
LINES 



(19) 



As pointed out in [T5HT7] . one expects the true 
JIMWLK result for higher point correlators of Wilson 
lines to significantly deviate from the simple product of 
dipole correlators (called "naive large N c limit" in [17] ) 
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used in the phenomenological applications so far (531 E3] • 
It was conjectured in Ref. [17], based on numerical tests 
for particular coordinate configurations, that JIMWLK 
evolution for the higher point correlators of Wilson lines 
can be approximated by a "Gaussian approximation" . 
Here one constructs, as in the MV model, the Wilson 
lines from color charge densities whose probability distri- 
bution is assumed to be local in rapidity and Gaussian. 
This enables one to express all higher point functions in 
terms of the two point function. The equations relating 
different Wilson line correlators are the same in the Gaus- 
sian approximation of JIMWLK and the MV model, the 
difference between the two being the different functional 
form of the two point function (dipolc). A theoretical ex- 
planation for the success of the Gaussian approximation 
is given in Ref. [151 H5]. 

The cross section formula ([5| involves correlators of up 
to six Wilson lines. The general Gaussian approximation 



for the six point function is not known, so we will 
here follow the approach used in Ref. [17] for the six point 
function. This approach consists of expressing all the 
higher point functions as products of fundamental repre- 
sentation traces, and then assuming that the expectation 
value of a product of traces factorizes into a product of 
expectation values. For the product of two traces appear- 
ing on the r.h.s. of the BK equation this has been verified 
numerically to be a very good approximation |44j . This 
"factorized Gaussian" has the advantage of preserving 
all the "coincidence limits" , i.e. it correctly incorporates 
the constraints coming from the SU(iV c ) group defini- 
tion V (x^V^ (Yt) —> 1, when x*r — > yr- While for the 
four-point function S^ 3 ' it would be possible to use the 
full Gaussian approximation, we will, for consistency, use 
the same approach also for S^ 3 \ Thus the formulae for 
the Wilson line correlators used in Eq. ^ are 



N 

S^(b T ,b' T , XT , x ' T ) » -^—S(x T ,x^) 



^UPsO^T, t>T) X T> x t) 

S( 3 >(b T ,x T ,z^) 



N 



Q(b T ,b' T ,x T ,x T ) - — g 



5(x T ,x^) 2 -^ 



N 



N c - 1 



S'(b T ,x T )S'(x T ,Zy) - —-^S(h T ,z' T ) 



(23) 
(24) 
(25) 



where we use the exact Gaussian quadrupole here. 

Q(hx, h' T , xj-, x' T ) calculated in Ref. [23]; since its In the large N c limit the Gaussian approximation re- 
expression is rather cumbersome we will not repeat it duces to 



5( 4 »(b T ,b^x T ,x^) 



5(x T ,x^) 



S(b T - x T )S*(x^ - h' T ) 



F(b T ,x T ,x^,b' T ) 



./V c ->-oo 



F(b T ,x^,x T ,b' T ) 
S(b T ,b^)S(x T ,x^) 2 



(5(b T - x T )S(x^ - b^) - 5(b T - b' T )S(x! T - x T )) 



S( 3 )(b T ,x T ,z^) « S(b T ,x T )S(x T ,z^) 



N, 



with the auxiliary function 



F(b T ,x T ,x^,b^) _ lnS{b T ,x' T )-lnS(b T ,b' T ) 
F(b T ,x^,x T ,b^) ~~ lnS(b T ,x T ) -ln5(b T ,b^) 



In S(x T , K) ~ 1q S(x t , Xy) 
lnS(x' T ,b^) -lnS(x^,x T )' 



(26) 

(27) 

(28) 
(29) 

(30) 



The first term in Eq. ( 26 ) is the "elastic" term, which is times referred to as the "inelastic" one. Note that one can 
the only one kept in Ref. [22]. The second term is some- explicitly verify that the limiting behavior — » Sp PS 



p p9 = 2 QeV, p a T ss = 1 GeV, yi = y 2 = 3.4 



p% 1 9 = 2 GeV, = 1 GeV, Vl = y 2 = 3.4 




3 4 
Atp [rad] 



uon parton level azimuthal correlation 
near forward RHIC kinematics, 



FIG. 3: The quark-; 
in dAu collisions, Eq. (17 1, 
at transverse momenta p^ ls = 2 GeV, s = 1 GeV and y = 
3.4. Shown are the "naive large N c " approximation with only 
the "elastic" contribution used in Ref. [23] and our Gaussian 
approximation of Eq. ( 23 1 and its large N c limit Eq. ( 26 1 



Note that a fixed A^-independent pedestal of 0.0002 GeV -4 
has been added to the "elastic" approximation for purposes 
of visualization. The DPS contribution is not included here. 



holds both for the factorized Gaussian approximation 
( 23 1 and its large N c limit ( 26 1 . The nonzero contri- 



bution in the DPS limit comes only from the "inelastic" 
part of which was neglected in Ref. [25] ; thus the log- 
arithmic divergence in this limit did not appear in that 
calculation. 
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VI. RESULTS 

We begin this section by showing results at the parton 
level, demonstrating the systematics as a function of dif- 
ferent trigger and associate momenta, different probe and 
target species. We then include the fragmentation func- 
tions for a more realistic comparison with both STAR 
and PHENIX dipion correlations. 

Figure [3] shows the effect of the "inelastic" contribution 
neglected in Ref. [35] to the dihadron correlation. Keep- 
ing only the "elastic" term, as done in Ref. [35] , results in 
an away side peak that is smaller by a factor of ~2 than 
the factorized Gaussian approximation, while the effect 
on the width of the peak is much smaller. Our conclu- 
sion from this plot is that including both contributions is 
essential for a meaningful comparison with experimental 
data. We also notice that the finite N c corrections to the 
Gaussian approximation can be neglected, and from now 
on we only use the large N c version of the Gaussian ap- 
proximation. The initial saturation scale for the nucleus 
is taken as Qf = 0.72 GeV 2 . 

Figure [4] illustrates the difference between deuteron- 



FIG. 4: The quark-gluon parton level azimuthal correla- 
tion at forward RHIC kinematics in proton-proton, proton- 
nucleus and deuteron-nucleus collisions. Top: normalized by 
the pedestal contribution, bottom: pedestal contribution sub- 
tracted (in which case the pAu and dAu results are identical 
by construction). Shown is only the large Ac result. 



nucleus, proton-nucleus and proton-proton collisions in 
the ratio of the peak to the A</?-independent background. 
Plotted on top is the total dihadron cross section, includ- 



ing the DPS contribution as discussed in Sec. IV divided 
by the minimum of the Ac/3-distribution. On the bot- 
tom this pedestal is subtracted. The initial saturation 
scales are taken as Q s q = 0.2 GeV for the proton and 
Q 2 = 0.72 GeV 2 for a nucleus. Comparing the pp and 
pA/dA results one sees clearly the depletion of the peak 
relative to pp collisions when the target is a large nu- 
cleus. This is the crucial generic feature whose observa- 
tion in the experimental data supports the concept of a 
depletion caused by intrinsic gluon transverse momenta 
in the target given by a saturation scale. Since we are ne- 
glecting any differences in the large x quark distribution 
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FIG. 5: The pedestal-subtracted quark-gluon parton level 
azimuthal correlation at forward RHIC kinematics, at dif- 
ferent values of the trigger transverse momenta p^ lg = 
1.5, 2, 2.5 GeV and y — 2.4. Shown is only the large N c result. 
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FIG. 6: The total quark-gluon parton level dihadron pro- 
duction yield at forward RHIC kinematics at different values 
of the trigger and associated particle transverse momenta at 
y = 3.2. The DPS contribution is included and the result is 
normalized by the pedestal yield. Shown is only the large N c 
result. 



between a deuteron and a proton, there is no difference 
in the correlated peak between pAu and dAu collisions. 
Due to the enhanced DPS contribution in deuteron col- 
lisions (where the possibility of taking one large x quark 
from the proton and the other from the neutron makes 
it much easier to observe a double scattering event in 
very forward kinematics), the DPS contribution is much 
larger in dAu collisions. This leads to a smaller ratio of 
the correlated peak to the pedestal, purely because of the 



FIG. 7: The total quark-gluon parton level dihadron produc- 
tion yield at forward RHIC kinematics, for trigger transverse 
momentum pl^ lg = 2 GeV, and different values of the asso- 
ciate transverse momentum pf? s = 0.5, 1,1.5 GeV and y — 3.4. 
Shown is only the large N c result. The DPS contribution is 
included and the result is normalized by the pedestal yield. 
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FIG. 8: The it azimuthal correlation compared to the 
PHENIX [TI] dAu result for two different p T rig bins. The 
Aip-independent pedestal in the plot is adjusted to fit the ex- 
perimental data, see Table|T]for the calculated estimates. The 
initial saturation scales are Q b q = 1.51 GeV 2 (solid line) and 
Q 2 = 0.72 GeV 2 (dashed line). 



increase in the denominator. The effect is similar to what 
is seen in the STAR forward neutron tagged (effectively 
pAu) data [45] and discussed in Ref. [40] . 

Figure [5] shows the dependence on the parton level 
dAu result on the trigger transverse momentum. Shown 
is the pedestal-subtracted yield per trigger for associate 
parton transverse momentum pt = 1.5 GeV and trig- 
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P, 3 < y h 2/2 < 3.8, 0.5 GeV < p a T ss < 0.75 GeV 



d + Au, 2.4 < y u y 2 < 4, 1 GeV < pf 3 < pp 9 , p% ia > 2 GeV 
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FIG. 9: The n azimuthal correlation compared to the 
PHENIX [11] pp result, for the same transverse momenta as in 
Fig. [8] The Aip-independent pedestal in the plot is adjusted 
to fit the experimental data, see Table [I] for the calculated 
estimates. The initial condition is the MV model with initial 
saturation scale Q 2 = 0.2 GeV 2 . Data points at larger p!^ ls 
are sifted by —0.03 GeV -1 to avoid overlap. 
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FIG. 10: Integrated yield under the away side peak in central 
dAu collision divided by the corresponding yield in pp com- 
pared to the forward rapidity part of the PHENIX data [11] . 



ger transverse momenta px = 1.5, 2.0, 2.5 GeV. These 
should be compared to the typical nuclear saturation 
scale Q s ss 1.5 GeV at these rapidities (defining Q s via 
the saturation criterion S(r — 1/Q S ) = 1 — e -1 / 4 ). One 
again sees clearly the expected effect of the nuclear sat- 
uration scale: the peak is small for small trigger trans- 
verse momentum, but starts to reappear when the trigger 
momentum is increased. Note that the rapidity in this 
plot is slightly less forward than elsewhere. This is to 



FIG. 11: The if azimuthal correlation compared to the 
preliminary STAR [10] result. The initial saturation scales are 
Q 2 = 1.51 GeV 2 (solid line) and Q 2 = 0.72 GeV 2 (dashed 
line). 



avoid hitting the kinematical limit x g + x q < 1 at large 
x in the deuteron. At very forward rapidites and high 
enough transverse momenta the kinematical limit causes 
the height of the peak in the yield per trigger to de- 
crease when the trigger pt is increased. This is caused 
by the fact that the correlated quark-gluon pair in the 
numerator of the ratio requires a larger total longitudi- 
nal momentum (and is thus suppressed at x — > 1) than 
the single quark (at the same y and px as the trigger) 
in the single inclusive cross section in the denominator. 
Because of this effect the dependence of the peak height 
on the trigger pj< at forward RHIC kinematics is actually 
smaller than the ratio of the momentum to the saturation 
scale would suggest. 

Figures [6] and [7] show the effect on the yield per trigger 
of varying also the associate pr- Since the two parti- 
cle yield (just like the single inclusive associate particle 
yield) is a steeply falling function of the associate pr, the 
yield per trigger itself falls as a function of pf? s . To show 
what is actually happening to the peak with respect to 
the background we again include the DPS contribution 
to the pedestal and divide the correlation by the Aip- 
independent part. As expected, the away-side peak is 
strongly suppressed for smaller trigger px- 

Figures [8] and [9] compare our calculation of the away 
side correlation to the PHENIX [11] data for different 
momentum bins. Figure [9] shows the proton-proton and 
|8]the central deuteron-gold result. For the latter we show 
calculations using two initial values for the saturation 
scale: Q 2 o = 1-51 GeV 2 , which would be the natural 
estimate for these central collisions, and a smaller value 
Ql = 0.72 GeV 2 , which would be preferred by some of 
the single inclusive data. 

Due to the significant uncertainty in the single inclu- 
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Data 


Pt range 


pedestal 


exp. 


PHENIX pp 


1.1 GeV <p T rig < 1.6 GeV 


0.04 


0.095 


PHENIX pp 


1.6 GeV < P ¥ g < 2.0 GeV 


0.02 


0.079 


PHENIX dAu 


1.1 GeV <p t T rig < 1.6 GeV 


0.10 


0.176 


PHENIX dAu 


1.6 GeV < p T rig < 2.0 GeV 


0.08 


0.163 


STAR dAu 


2 GeV < p^ is , 1 GeV < p^ s < p^ is 


0.02 


0.0145 



TABLE I: Calculated estimates for the pedestal height compared to the experimental values. In PHENIX results the units are 
GeV -1 . The dAu values are for central collisions, and Qi is taken as Q 2 — 1.51 GeV 2 . 



sive baseline spectrum, the Ay-independent pedestal val- 
ues are rough estimates that cannot be directly com- 
pared to the data at this point. For plotting purposes 
the pedestal values have been adjusted to the data in the 
plot; the results of the calculation compared to the ex- 
perimental pedestal are given in Table]!] The normaliza- 
tion uncertainty also affects the height of the away side 
peak, even though we would expect that at least part 
of this theoretical uncertainty cancels in the yield-per- 
trigger ratio. Taking into account these uncertainties, 
we obtain relatively good description of the back-to-back 
peak, underestimating the peak height especially with 
Q 2 = 1.51 GeV 2 for the nucleus and in the case of the 
proton-proton scattering. 

PHENIX has also published data for larger trigger par- 
ticle momentum 2 < p£ lg < 5 GeV. We can not de- 
scribe the observed relative increase of the peak height 
compared to lower values of the trigger momenta in this 
kinematics. This phenomenon can be seen in parton level 
results (see Fig. [5]) , but in PHENIX kinematics we are so 
close to the kinematical boundary that in our calculation 
the deuteron parton distribution function suppresses the 
peak. 

Figure [10] compares our result for the nuclear mod- 
ification factor JdAu of the area under the peak to 
the PHENIX experimental result [IT]. The exper- 
imental data are given as a function of Xf rag = 

((^ rig )) e - + (24 ss ))e- /y^jv. We compute 

this quantity by calculating the ratios of the areas under 
the pp and dAu away side peaks at various kinematical 
configurations at rapidities 3 < y 1,1/2 < 3.8 and at var- 
ious trigger (1 < p^ lg < 2 GeV) and associate particle 
(0.5 < p^ s < 1-5 GeV) transverse momenta. We then 
obtain JdAu by averaging the results in every a;f rag bin. 
This is not exactly the same method that is used in ex- 
periments, but we have tested that the result does not 
depend much on the kinematical limits. Our calculation 
agrees with the PHENIX result within errors, and we 
clearly see that the suppression increases when Xf rag de- 
creases, which is expected as the saturation effects should 
become more visible at smaller Xf rag . The peak height is 
underestimated in both pp and dAu collisions probably 
due to the uncertainties in the single inclusive baseline 
normalization, and we expect that this error partly can- 
cels in JdAu. 



STAR [TU] result. While the errors in the data are rather 
large, let us point out a few things. Our calculation seems 
to agree reasonably well with the data, although perhaps 
underestimating the hight of the peak. It is clear from 
Fig. [3] that neglecting the "inelastic" part of the cross 
section (as done by Albacete and Marquet [35]) would 
make the agreement worse. 2 Note also that we are as- 
suming collinear factorization of partons into hadrons, 
which is very crude at these small transverse momenta. 
Any /c<r-smearing from fragmentation would broaden the 
away-side peak. 



VII. CONCLUSIONS 

We have in this paper performed a calculation of nu- 
clear modifications of forward dihadron correlations in 
the CGC framework. We use a running coupling BK 
evolution for the dipole and a factorized Gaussian ap- 
proximation for the higher point functions of Wilson lines 
to describe the target. Our calculation is the first one in 
the literature to evaluate the the full expression of the 
dihadron cross section, keeping both the "elastic" and 
"inelastic" terms (unlike Ref. [35]) and not restricted to 
the high-pT "correlation limit" (unlike Ref. [25J). We find 
that including the "inelastic" term in has a significant ef- 
fect on the correlation, enhancing the peak by a factor 
~2. We also show that in the appropriate kinematics the 
dihadron cross section reduces to Ay-independent dou- 
ble parton scattering, which in the earlier literature has 
been considered as a completely separate contribution. 

We show explicitly how the large saturation scale in 
a nucleus leads to a smoothing of the away-side peak 
in dihadron correlations for transverse momenta of the 
order of the saturation scale. We obtain a reasonable 
ordcr-of-magnitude estimate also for the Ay- independent 
pedestal part of the correlation. Taking into account 
the lack of a parametrization of the dipole cross section 
that would simultaneously reproduce both mid- and for- 
ward rapidity single inclusive hadron spectra from all the 
RHIC experiments (including the STAR and PHENIX 



Finally Fig. 11 compares our calculation to the 



2 When comparing to the result in Ref. 1351 note also that there is 
a significant numerical error in the calculation of Ref. [351. which 
we thank C. Marquet for indicating to us. 



12 




Acknowledgements 

We thank K. J. Eskola, I. Helenius, R. Paatelainen, 
B. Schenke, M. Strikman and R. Venugopalan for discus- 
sions and J. Albacete and C. Marquet for helpful com- 
parisons with their results. H.M. is supported by the 
Graduate School of Particle and Nuclear Physics. This 
work has been supported by the Academy of Finland, 
projects 141555 and 133005, and by computing resources 
from CSC - IT Center for Science in Espoo, Finland. 



FIG. 12: Illustration of the kinematics in Eq. |A5l 



Appendix A: S (4) and S <2) in the DPS limit 



forward ir° data sets for which the correlation measure- 
ments are performed) it is difficult to consistently per- 
form a more accurate calculation of the pedestal contri- 
bution at this point. We thus leave a more thorough ex- 
ploration of the single inclusive spectra for future work. 
Results from the future LHC proton-nucleus run, with 
the significantly larger kinematical coverage available, 
should significantly clarify these uncertainties. 



Let us now show that only the four point function 
gives an infrared divergent contribution to the the double 
inclusive cross section. We then show how this contribu- 
tion can be identified as the double parton scattering con- 
tribution appearing in the literature, although calculated 
in the "hybrid formalism" and not collinearly factorized 
perturbation theory as usual. 

In the DPS limit and for massless quarks the double 
inclusive cross section Eq. ^ reduces to 



dy 9 dy 9 d 2 q T d 2 k T dps 



xq h (x)z(l - z)a 



d 2 x T d 2 x^ d 2 b T d 2 b' .(x^- 



Pn 



(2tt) 2 (2tt) 2 (2tt) 2 (2tt) 2 
qg (z) \s(b T - W T )S(x T - x^O H 



XT) e iqr-(by-b T )g 7r 2 



ut| 2 K p 



S(z(x T -x' T ) + (l-2)(b T -b' T )) , (Al) 



where we recall that S denotes the adjoint representation 
two point function and that ut = xt — by; u' T = x' T — 
h' T . We have denoted the splitting function for radiating 
a gluon from a quark by P q ^ qg (z) — Cf(1 + (1 — z) 2 )/ z. 
Using the momentum representation 

^^/^ 5' (A2) 



shifting the integration variables xj- and by by x^ and 
h' T respectively, then integrating over x' T and h' T we get 



fo-hA-yqgX 



dy 9 dy 9 d 2 q T d 2 k T dps 



xq h (x)z(l-z)a s S ±7 ^ f d 2 x T d 2 b T f n L 2 e-^~ nT ^ T e-^ T+n ^ b ' 
{2n) b J n T 2 



P q .+ qg (z) S(b T )S(x T ) + S(zx T + (1 - z)b T ) , (A3) 



Let us first look at the second term, which originates 
in S^ 2 K Changing the integration variable from by to 
zt = zxt + (1 — z)hx the z<r-integral gives a two-point 
function in momentum space and the x^-integral a 5- 



function which can be used to integrate over iif to get 
f \ /t ^ a 1 od / \ S(q T + k T ) 

xqh (x)z(i - «Mi^2PH«Wpr^r^- 

(A4) 
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This is a finite result, which is included as such in 



Eq. (17) 



The first term of Eq. ( A3 1 , on the other hand, becomes 

- — - — ^ jk. — = xq h (x)z(l - z)P q ^ qg (z) 

dy 9 dy 9 d 2 q T d 2 k T dps 



a s S± 



1 



(2tt) £ 



d 2 riT 
riT 2 



2S(q T + n T )S(k T - n T ). (A5) 



The interpretation of this result is illustrated in 12 the 
emission a gluon with transverse momentum is fol- 
lowed by the subsequent independent scattering of the 
quark and the gluon off the target. The integral over 
riT is logarithmically divergent in the infrared iiy — > 0. 
Regulating it at a parametrically confinement scale mo- 
mentum A we can neglect compared to kr, qr in the 
two point functions and approximate 



da hA ^ qgX 



dy 9 dy 9 d 2 q T d 2 k T dps 



= xq h (x)z(l - z)P q ^ qg (z) 



^InA 2 1 
2tt 



rS(q T )S(k T ). (A6) 



(27T) 



We find it instructive to compare this to the general form 
for the LO DGLAP evolution equation for the quark- 
gluon double parton distribution D qg (x q , x g , Q 2 ) (see e.g. 



Eq. (2.1) of Ref. [15]): 

Q 2 ) a s (Q 2 ) r 



ding 2 

q(Xg +Xg,Q 2 )- 



2tt 



evol. 



P, 



(A7) 



Here "evol." denotes DGLAP evolution terms corre- 
sponding to separate radiation from quarks and gluons, 
leading to a variation of the double parton distribution 
with the scale Q 2 . As we are not interested in very high 
transverse momentum scales we will neglect these and 
concentrate on the second term, which corresponds to 
the creation of a new correlated quark-gluon pair in the 
projectile due to radiation from the quark. This leads us 
to identify 



D qg (x q: x g ,Q ) 



a s (Q 2 ) 



2ir 



In A 2 q(x q 



1 



x 



-P, 



,1 a 



uq t Xg \ ~r ^ g 

With this identification wc write the contribution of 



q^qg 



x 



X r . 



(A8) 



Eq. (A6) as 



foh-A-yqgX 



dy g dy s d 2 q T d 2 k T dps 



z(l—z)x 2 D qg (zx, (1— z)x, A 2 ) 



S± 
(2tt)< 



5(q T )5(k T ). (A9) 



In the uncorrelated approximation where the double par- 
ton distribution is approximated by 

D qg (x q , x g , Q 2 ) x q q{x q , Q 2 )x g g(x g , Q 2 ), (A10) 



the result of Eq. ( A9 1 can be identified as the product of 
the single inclusive cross sections |l]) and ^ as 



j HA— >qqX 

dcr DPS 



1 da hA ^ qX da hA ^ gX 
dy q dy g d 2 q T d 2 k T S± dy q d 2 q T dy 9 d 2 k T ' 



(All) 



The prefactor, which must have the dimensions of an 
inverse area, is commonly denoted as l/cr ff- In our cal- 
culation, which assumes a dilute probe hitting a large, 
uniform target, this misleadingly named "effective cross 
section" is the same as the transverse size of the target 
S±. 



Appendix B: Impact parameter profiles 

Let us start (see e.g. Ref. [42]) by considering the gen- 
eral case of two colliding projectiles 1 and 2, of partons 
distributed with an impact parameter profile that we will 
assume is factorized into a distribution Fx 2 (br) normal- 
ized as 



J d 2 b T F h2 {b T ) = 1. 



(Bl) 



We define the Fourier transform as 

F(A T ) = J d 2 b T e lbT AT F(A T ). (B2) 
A single inclusive cross section is simply 

<t s = J d 2 b T1 d 2 b T2 F 1 (b T1 )F 2 (b T2 )^ 

= J d 2 b T d 2 b T iFi(b T i)F 2 (b T -b T i)cr lJ ', (B3) 

where a 11 is the parton level cross section multiplied by 
the appropriate parton distribution (which we assumed 
factorized from the impact parameter profile) and b^ the 
impact parameter of the collision. For a double scattering 
the cross section can be written as 

a D = J d 2 b T d 2 b T1 d 2 b' T xFx(b T i)Fx(b' T x) 
F 2 (b T ~ b T i)F 2 (b T - b' T1 )a ij a kl 
= J |^l^i( A T)| 2 |F 2 (A T )|W fci . (B4) 
The term "effective cross section" refers to the factor 

— = / U%I^i(At)| 2 |^(A t )| 2 , (B5) 

and has experimentally been measured as <7 ff f=a 15 mb 
in proton-antiproton collisions [171 HS]- The value 15 mb 
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is surprisingly small (i.e. the cross section surprisingly 
large) compared to other estimates of the proton im- 
pact parameter profile, which could lead one to ques- 
tion whether it includes some additional dynamical cor- 
relations that should properly be included in the dou- 
ble parton distribution itself. For example a Gaussian 
parametrization of the TOTEM [49J B = 23.6 GeV" 2 
elastic cross section da / dt ~ e 23 6t / GcV would corre- 
spond to F(A T ) = e -5-9^T 2 / Gov 2 and CTcff > 100 mb 
A perhaps more relevant data point for hard small x 
processes would be diffractive vector meson production 
in DIS [50j [51], where we can parametrize F(At) — 
e -i3 D A T 2 /2 with Bd ~ 4 GeV" 2 , leading to a eS « 39 mb. 
The dipole parametrization estimate |52j for the proton 
two gluon form factor, F 2g (A T ) — 1 / {A T 2 / m 2 g ) + l) 2 
with m 2 = 1.1 GeV 2 , similarly leads to a large value for 

Ccff- 

It was, however, noticed in Ref. |40j that the experi- 
mental result 15 mb could be interpreted as the large x 
valence-like partons being completely localized in the 
transverse plane. In this case the dilute probe form factor 
would be Fi(At) = 1 and 



Equation ( B6 ) leads to an intuitive picture of the double 



— = / ^%I^(A T )| 2 = / d 2 b T F 2 (b T ) 2 



inclusive cross section as an integral over the transverse 
plane of the target that is directly generalizable to nu- 
clei. For a nuclear target with ^(bx) ~ Ta^t) this 
would lead to a similar impact parameter dependence 
as the contributions (b) and (c) in Ref. [40]. Here the 
target is proportional to Ta(\>t) only in the dilute limit 
and assuming that the target density is ~ Ta would be 
inconsistent with a nuclear modification factor R p a sig- 
nificantly different from one in single inclusive scattering. 
However, motivated by the observation in Ref. [ID] that 
the proton-proton result of a e s ~ 15 mb is consistent 
with the large x projectile being pointlike we shall keep 
the interpretation (B6| for <r c g. Thus the only impact 



parameter dependence in our calculation is an integral 
over the transverse profile of the target, which disappears 
when calculating yields in stead of cross sections, as in 
Eqs. Q,©, and @. 



(B6) 
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